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Abstract 

We discuss experimental signatures capable of nearly immediate study that 
would discern /constrain new physics manifested via enhanced gluonic penguin 
decays of the b. 
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I. INTRODUCTION 



Despite the large B mass, it appears that a systematic expansion in a s and l/m^ fails 
by about 15% PJ to predict the semileptonic branching ratio 



Br s, = i^-^J {) 
1 B 

This discrepancy presumably requires enhancement of the nonleptonic decay width and most 
likely will be resolved within the framework of the Standard Model. For example, it could 
be that local "quark-hadron duality" is not as good at the B scale as originally hoped. 
In particular, one recent suggestion is that non-perturbative effects could significantly 
enhance the b — > ccs rate (with relatively slow charm quarks in the final state) beyond the 
naive T(b — > ccs)/Y{b — > cud) 1/3 phase space prediction. 

Alternatively it has been pointed out that an enhanced b —* sg rate via penguin ampli- 
tudes could help resolve the discrepancy. 0] These will contribute only to the nonleptonic 
decay width, r . These decays do not lead to charm production (and hence the reason that 
their contribution could be added incoherently). Therefore such enhanced penguin decays 
would not generate an unwanted charm excess [0. 

A branching of T(b — > sg)/^ 1 « 20% is needed to explain completely the semileptonic 
deficiency. To leading order in a s , the gluonic penguin decay rate is given by 

nb ^ sg) A(Y^MM\\ (2) 

7T \ Vbc J 

To = m foG ! FHc/1927r 3 is the lowest order, perturbative result for the semileptonic decay 
b — ► cez/ e , (excluding phase-space modifications due to the finite mass of charm quark), and 
Cm(i^) is the coefficent of the magnetic penguin operator entering the effective Hamiltonian 

(Heff) 



Of = ^ m b sr a G%a^ (1 + 75 ) b, (3) 



evaluated at a mass scale \i appropiate to the decay at hand. Working within the standard 
model and taking fi = m&/2, one obtains, using the evolution equations for Cm found in 



Ref. (wherein the recent controversy ||,^| concerning the apparent scheme dependence of 
this evolution has been resolved) that 

Br(b ^sg)^2x 10~ 3 . (4) 

We have taken V ts = .045, the lifetime t b = 1.3 picoseconds J7| and we have equated the 
bottom quark mass with that of the B meson, i.e., rrib = Mb- Using m;, = Mb — .5 yields a 
rate 60 % of the above. 

A significant enhancement above standard model predictions is thus required for gluonic 
penguin decays to play a role in the semileptonic decay rate issue. Admittedly this scenario 
may seem unlikely in view of the recent measurement by CLEO || of the inclusive photo- 
penguin decay rate, b — > 37, found with a branching rate of only Br(b — > 57) m 2.2 x 10~ 4 , 
and in general agreement with Standard Model predictions ||. Nonetheless, a preferred 
enhancement of the gluonic penguin by some appropriate SUSY extension of the Standard 
Model relative to its electromagnetic analogue is conceivable. Even an enhanced branching 
ratio Br(b — > sg) of just a few percent could be part of a "cocktail" [0 solution for the 
semileptonic problem. 

In the following we argue that inclusive measurements of the kaon spectrum (presumably 
Kg for experimental feasibility) near the region of u = Pr / PRmax ~^ 1 would be sensitive to 
b — » sg rates in the range 5% — 20%. We also discuss exclusive decay modes, in particular 
B — > Kit, first indications of which were possibly seen last year at CLEO |1(J. Here too 



enhancements would naively be likely, although as we will see, bound state effects could 
complicate such expectations. 



II. INCLUSIVE K SPECTRUM 

For a B meson at rest, the inclusive kaon momentum spectrum for kaons emerging from 
the cascade 
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B -> D + X 

^ K + X' (5) 
can be readily calculated from the available data. The inclusive momentum distributions 



for B — > D(P) + X [11] and D — ► + X' [12| have been experimentally measured. A 



good overall fit to these distributions can be given by the rather simple parametrizations 

Pd/b{x) = 6s(l - x) 
6 

Pk/d{v) = J7-2 ^(V - Vmin){l - y) (6) 

where 

= P / Pmax V = Q/ Qmax Vmin ~ -14 (7) 

and all momentum are defined in the rest frame of the decaying heavy meson. Note that 
the crucial regions of interest, x, y — *■ 1 are particularly well fit by Eq. (|6|). For simplicity 
we have also normalized each distribution to one. 

The energy Ex of the kaon in the B decay frame is obtained from the energy and 
momentum in the D decay frame via a boost transformation: 

M D E K = P°q° + \P\\q\z (8) 

where z is the cosine of the angle in the rest frame of the D meson between the decay 
direction of the kaon and the boost direction of the D. Defining 

u = E K /E Kmax , (9) 



and 



a 



= Mb 2M b M d 
P max M 2 _ M 2 

_ = 2M D M K 

p maxqma , _ {Ml - Mp{Ml - Mp 

M D E Kmax 2M 2 D {M 2 B + M 2 D ) 1 J 



we have then that 



u = a 



V x 2 + ot 2 \Jy 2 + f3 2 + xyz 



(11) 



Using the fact that the D is spinless so that the angular distribution z is uniform, one 
obtains that the inclusive energy distribution Pk/b{u), of the cascade kaons is given simply 
by 



Pk/b(u) 



dxp D / B ( 



x 



1 ,1 

dyp K /D{y) / dz8(a 



f-J rrt i n 



V x 2 + a 2 J y 2 + (3 2 + xyz 



u). (12) 



Integrating over the 5-function using the z integral, Eq. ([TJ) becomes 

Pk/d{v) i 



Pk/b( 



u 



1 dx pm^lf 1 dy > 



6(a Vx 2 + a 2 ^y 2 + (3 2 + xy — u) x 



Ofw — a 



\/x 2 + a 2 \Jy 2 + /9 2 — 



(13) 



Figure (1) displays the resulting distribution Pk/b{u) which for conformity with the rest of 
the literature, |TT| , |12"|1 , we have plotted not as a function of the energy variable u, but in 
terms of the kaon's momentum, 

U = P K /PKmax- (14) 

The two remaining integrals in Eq. (|T3|) were performed numerically and the entire distri- 
bution was again normalized to integrate to one. 



We note that at the present only preliminary data |18[ on the Kg inclusive spectrum at 
the T As is available. Nevertheless this does conform to the expected overall cascade form 
of Figure (1). Our main interest here however is in the endpoint region, u ~ u — > 1 for 
which much more precise data will be required and should be available from CLEO II. In 
this endpoint region one can directly show that 



\\-mp K/B {u) oc 

u— >1 



l-U)\ 



(15) 



For u = 1, both momentum variables x, y must also approach their upper limit. The 



explicit distributions Pd/b(x) and pK/B(y) thus yield two of the four powers in Eq. (|l5l). 



The remaining two powers arise from phase space constraints generated by the functions 
which force x,y > 1 — 77 if it = 1 — i] and rj — > 0. The importance of this result is that 
Pk/b{u) is severely suppressed near u — 1 compared to the much harder function for an s 



quark jet to fragment into a kaon. From general counting rules fOUII^I one can show from 



perturbative QCD that the fragmentation function for ans-> K +, ° behaves as 

lim D K+ , 0/s {u) ~ (1 - uf. (16) 



u-»l 



Combining Eqs. (|l3f) and (|T6|) , the complete distribution, p^* B (w), of kaons from 5 
decays, is given by 

P*/b(«) = (1 - c)Pk/b{u) + e/^+.o/^w) (17) 

where e is the b ^ sg branching fraction. From the previous discussion, one knows that 
for any finite e, at some sufficiently large u — > 1, the s quark fragmentation function must 
dominate. The issue now is to better quantify this result with a reasonable parametrization 
for D K +,o /s (u). 

The new CLEO B — > X s + 7 data [^TJ could in principle yield s ^ K information in a 



setting which appears to be kinematically similar to that in b — > sg. Nearly 30% of all the 
B — > X s + 7 inclusive decays occur through the K* resonance. If we simply focus on the 
kaons produced through the decay of the K*, one finds that the kaon spectrum is very hard 
as u — > 1, with a typical value of u = .7 (replacing M D by in Eq. (p]), taking x, y ~ 1, 
and noting that z ~ for the p-wave decay K * — > Kit) . Such a contribution would show up 
quite dramatically in P^/ b (m) for a penguin decay rate of order 10% and indeed would allow 
ready determination of e in Eq.([T7|) to values as small as e = .02 by integrating over the 
total kaon yield above u = .7. In an ideal case of infinite B mass the s quark from b — > sg 
or b —>■ S7 would have identical fragmentation independent of the identity of the recoiling 
system (gluon jet or single photon). However the use of this inclusive data for the actual 
B mass, Mb = 5.3GeV, is somewhat dubious. For b — > 57, the relevant total energy of the 
hadronizing system of the s plus spectator q is rather low, being roughly only 
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W Xs = ^2E s a m q " « (5.3GeV x .3GeV) 1/2 = 1.3GeV } (18) 

as indeed is manifest by the fact that the K* and only a few other kaon resonances domi- 
nate the data. Barring precocious scaling, extraction of accurate information on inclusive 
functions such as Dk+,o/ s {u) in such kinematics is questionable. 

The kinematics in the case of b — > sg is, on the other hand, different. Assuming the s 
quark to recoil against an oppositely moving gq = 3 source, the total energy available for 
hadronization is now Mb = 5.3GeV and we expect the "leading" s — > K fragmentation to 
be similar to that in an e + e~ — > ss process at these energies. Indeed comparing even just 
the expected standard model rate, Eq. with the upper bound given last year by CLEO 
for B — > Kir ( < 2.6 x 10~ 5 ) indicates that the inclusive rate will not be dominated by just 
a few resonances. 

Unfortunately s — > K fragmentation contributions in e + e~ is charged suppressed by a 
factor q^/Y^Qi an d are hence difficult to extract from the known data on inclusive kaon 
production [TH|. Likewise, s quark jet production in deep inelastic v scattering is Cabibo 
suppressed by \V US ^ and hence not particularly useful. We are thus forced to work by 
analogy, guided by general symmetry principles. 

Starting first from £77(3) flavor symmetry, the inclusive 7r + spectrum spectrum from 
e + e~ annihilation has been recently [16| nicely fitted using a primary, direct u — > 7r + term, 



D p , and a secondary distribution term D s 

D p (z,t ) = l(l-z) 2 Ds ( Zjto ) = l<J—^ (19) 

o o z 

where to = ln(Ql/A 2 ) reflects the general fact that in QCD these distributions run with 
Q 2 . The motivation of the authors of |16| for departing from other, perhaps more common, 
forms in which D(z) is given jointly as one smooth function (e.g. as in ]T7| where D(z) = 
(1 — z) 2 /Az) is that the physics dictating the two end point regions is rather different. The 
logarithmic rise in total kaon number produced for z — > is driven by mesons produced 
out of secondary quarks formed in the fragmentation chain of the outgoing quark jet. The 
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region z — > 1 is dominated by mesons made out of the original quark in the jet. A separation 



of these two phenomena was crucial for the study in jl6] in which (as is the case here) the 
authors were particularly interested in the z — ► 1 regime of the fragmentation function in 
order to compare with competing processes. 

Ignoring completely secondary kaon production, and assuming for the kaon the same 
primary distribution D p (z,t) as in the pion, (evolved down from the fitted Qo = 29GeV 
data, to Q — Mb using the analytic expression provided in fl6|D, the result for p^/ B (w), Eq. 
(|17|), using a value for e = .2 is shown in Fig. (2). For comparison is included the result of 
Fig. (1) (i.e. e = 0), and as is appropriate, we have focussed only on the end-point region 
u — > 1. We see that for this "maximal" value of e, a significant difference has developed in 
the expected kaon distributions once u > .7. Accurate data at larger u allows one to probe 
significantly smaller e, so that for e = .05 comparable differences appear at u = .8 and at 
u = .9, a value as small as e = .02 could be discerned. 

Having demonstrated what should be an experimentally testable effect, it is likely that 
these estimates are however conservative. SU(3) flavor symmetry is in fact broken and it 
is likely that the s quark fragmentation function is harder than that of a u quark. As an 
extreme alternative, we show the results in Fig. (3) of using the Heavy Quark fragmentation 
function of Peterson et al. Jl9[ 



N 

D K+ ,o /s (z) = — — -j (20) 

z [1 - 1/z - e Q /(l - z)\ 

in which iV is determined by fixing the normalization to integrate to one and €q is qualita- 
tively m^/rriQ, the ratio of effective light to heavy quark masses. Note that the e — * oo limit 



smoothly matches up to the parametrization of Baier et al. |17 for a light-quark fragmen- 
tation function. In the case of charm a good fit to the data |2(| is found using ec = .15. 
Since in the case of the s quark the choice of e s is more ambiguous (if at all correct), Fig. 
(3) contains plots for a few possible e s values. As expected, the resulting distributions yield 
a significantly greater departure from the cascade scenario of kaon production than the one 
indicated by Fig. (2) where perfect SU(3) flavor symmetry was assumed. 
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III. EXCLUSIVE DECAYS 



The observation last year by the CLEO group of the two body exclusive decay, 
B — > K*^f was the first unambiguous experimental evidence of penguin processes. One 
likewise expects that rare two-body hadronic decays of the 5-meson to be a particularly 
useful means of measuring gluon mediated, h — > s transitions. For concreteness, we will 
focus on the mode B — > Kn which, as mentioned earlier, was likely seen last year at CLEO 
|T0| (the ambiguity involves insufficient experimental resolution to separate candidate im 
from K7r decay channels). 

For the two body hadronic decays, penguin processes compete with another rare decay, 
b — > u, which occurs at tree-level in the Standard Model but is proportional to the CKM 
matrix elements IV^V^ 2 . Judicious comparison with analogous processes proportional to 
|KifeKid| 2 allows one to infer the relative importance of the gluonic penguin contribution. 
In the case of B — > Kix the corresponding decay mode is B ^ ttti. (Hence the further 
importance that CLEO resolve these two modes, the sum of which were reported with a 
total branching rate of 2.4 ±'y ±.2 x 10~ 5 [|lOj) An observed ratio of branching rates much 
above (or below) the naive (Vus/Kczl 2 ~ 1/20 must be due to penguins. 



Using perturbative QCD methods recently seen [22] to give a good description of the two 



body hadronic decays of the B meson, we estimate in the standard model that 

Br(B ->• Kir) ss .5 x 1(T 5 , (21) 



not far from the CLEO data of last year [|K]]. Such a result clearly does not allow much room 



for enhancement. However, the importance of bound-state effects must be emphasized. We 
will therefore sketch how our estimate Eq. ( pl|) was obtained, in which some particularly 
simplifying approximations were used. A full analysis of the decay is the subject of a 
forthcoming work. 
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Since the two body exclusive decays of the B involve large mometum transfers, they are 



short distance events. A twist expansion in perturbative QCD suggests [23] that only the 
contribution from the lowest order Fock component expansions of the B and of the outgoing 
mesons are relevant. The decay rate of the B then involves a perturbatively calculable 
hard amplitude convoluted with a soft physics wave function, ip m , from each of the mesons. 
These wavef unctions, although as yet uncalculable from first principles, are universal for 
each meson, i. e. they factorize from the hard amplitude and hence are independent of the 
process involved. Thus as was employed in Ref. f22| , ideally one can phenomenologically 



parametrize these wavefunctions using a (few) measured cross-sections/decay rates. For 
simplicity, we use the factorization scheme advocated by Brodsky and Lepage [fj4j and take 



the momenta of the quarks as some fraction x of the total momentum of the parent meson 
weighted by a soft physics distribution amplitude (f>(x) (4>(x) being then simply the quark's 
wavefunction tp{x,k±) integrated over transverse momentum, k±). 

An important ingredient in the perturbative QCD approach to the two body exclusive 
decays of the B is that the decay amplitude can acquire an imaginary part because some 
heavy quark propagators in various Feynman graphs can go on-shell in the integration over 



the mesonic distribution amplitudes. Here as elsewhere P5|] , |2jJ, picking up such poles is 
legitimate in pQCD as they are not pinched singularites and hence by the Coleman-Norton 
theorem f27f , are not associated with a long distance event. These poles arise (in part) 
because of the factorization scheme we employ and because of what one believes to be the 
correct relation between quark and meson masses: 

M B = m b + A, (22) 



where A ~ 500MeV. Since it has been found in practice j22J that these imaginary parts tend 
to dominate the decay amplitude when they occur, we will here focus on them to obtain our 
estimate Eq. ( |2T]) with the more complete study to be presented elsewhere J23] . 



For our present purposes then, the most relevant graph contributing to the decay rate 
B — > Ktc is shown in Figure (4). The square indicates a gluonic penguin, the operator 
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structure of which will be presently discussed. The cross indicates the heavy-quark propa- 
gator that can go on-shell. We note that graphs at this order in a s other than that shown 
have imaginary parts, but they are suppressed by factors of e = A /Mb- Unlike the case 
of the inclusive cross-section considered in the previous section, the restriction to explicit 
hadronic modes in the final state means that the virtuality of the quark and gluon legs en- 
tering the penguin decay b — > sg cannot be discounted. Hence more than merely the simple 
chromo-magnetic penguin operator Off, Eq. @, contributes and in particular, it is found 
that chromo-electric penguins 



Of = if^V" (1 - 75) b (DpGrY , (23) 



are highly relevant (indeed they in fact dominate our estimate for the B — > Kir decay 
rate). Including then both operators Off and Of, with Wilson coefficents Cm(/i) and 
C#(//) respectively, the contribution (I) of the diagram in Figure (1) to the decay amplitude 
is given by the expression 

(/) = —c M Qi) I dx Mx l {1 ~ x) f dyMv)y 

eb J x — 2es — vt] J 

AA r fa(x)(l -x)(l + x- 2e B ) r 

+—C E (fi) / dx : \ dyfa(y)y(l - y) 24 

ts J x — 2e B — it] J 



where A is given by 



A = -a 2 J B fkUG F U bt U; s , (25) 



and a peaking approximation has been used for (f> B {z 

1 

2\?3 



B ( z ) = J-f B 5(z-e). (26) 



fs is the decay constant of the B. 

Using the distribution amplitude of Chernyak and Zhitnitsky |28[] for the pion and for 
simplicity, the asymptotic distribution amplitude |29j for the kaon, 
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4> v {x) = 5(1 - 2xf 

My) = i (27) 

the imaginary piece of Eq. ( f2"4|) becomes: 

j m (j) = (-C^/i) + 4C M (//))(1 - 2e)(l - 4e) 2 . (28) 

We use Im(I) to obtain our estimate of the decay rate 



16tcM b 



As parameters we use 



V t8 = .045 

/b = V2f n 
.5 

£ ~ Mb 

A QCD = -2GeV. (30) 

For the scale \i we take the virtuality of the softest gluon exchanged, /z 2 ~ .5(Ge^) 2 . For the 
evolution of Ce(/j), only Cella et al. of Ref. have calculated the anomalous dimension 
mixing matrix relevant for the chromo-electric penguins, and thus we use their results. 
This might be thought imprudent, considering the controversy that was associated with the 
chromomagnetic penguins [f§-|6|]. We note however that in the case of the chromomagnetic 
penguins, the final results of Ciuchini et al. M differ only slightly from that of Cella et al. 
for those operators mutually calculated, and hence one might expect that any errors due to 
scheme dependence would be kept at a minimum. Such is our hopes in the present work, 
although we acknowledge and stress the importance that these expectations be confirmed. 
With these inputs we obtain that 

Br\B -> Kn) w .3 x 10~ 5 , (31) 

where the superscript I reminds the reader this estimate is based only on picking up the 
(leading) imaginary part of the decay amplitude. Roughly assuming that the real part 
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(which adds incoherently but involves many more graphs and also in principle more off- 
shell operators) is of the same order, we obtain Eq. (pip . 



We note that similar such bound-state effects have been previously considered [30j in the 
case of B — > K*^ decays, however there the physical on-shell photon eliminates the elec- 
tromagnetic version of Eq. Only bound-state effects involving the off-shell character 
of the quark propagators were thus needed to be considered. Although large in amplitude, 
these amusingly were found, due to accidents of phase, to have minimal effect on the decay 
rate assuming standard model parameters of H e ff. The purely hadronic two body decay 
B — > Kit does not however enjoy such simplifying features. Chromoelectric penguins are 
found to be highly relevant because the restriction to a Kir final state allows for the gluons to 
be significantly off-shell and because the Wilson coefficient of the chromo-electric operator 
is appreciablely larger than that of the chromo-magnetic operator (at either the My/ scale 
or when ran to lower scales appropriate for the decay of the B). 

The complications of bound-state effects and in particular the introduction of additional 
operators thus allows, at least in principle, significant differences between the inclusive rates 
discussed in section (II) and that of any particular decay mode as discussed here. One ten- 
able although perhaps contrived scenario is that only Og is significantly enhanced. Although 
somewhat bizarre, especially since for the gluon penguins it is at a fundamental level the 
same Feynman graphs that determine the Wilson coefficients of both the chromomagnetic as 
well as chromoelectric operators, such a scenario cannot nevertheless be precluded. Indeed 
the fact that QCD corrections play an important role and produce significant changes in 
both the absolute magnitudes and relative sizes of the various coefficients, means that such 
a scenario might yet be feasible. However one should recall that this preferential treatment 
would also have to extend into the electromagnetic sector where the data [p|,pT| also does 
not allow significant departure from standard model predictions. 
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IV. CONCLUSIONS 



It is well known that deviation of b — > 57 and b — > sg amplitudes from standard model 
estimates could directly indicate new physics. The study of potential b — > sg enhancement 
is further motivated by its possible contribution to resolving the B semileptonic problem. 
We have here discussed the possibility of detecting such enhancement and estimating the 
overall strength of b — > sg transitions by focusing on the resulting inclusive kaon distribution 
in the region near Pk = p™ ax and alternatively looking at the extreme case of Pk = P K W,X 
corresponding to the (penguin generated) B — > Kit exclusive final state. 

Clearly future studies could make use of the much richer topologies in an effort to have a 
more sensitive extraction of b —>■ sg. [] Remarkably however, mere use of the inclusive kaon 
distribution will be sufficiently sensitive to an enhanced b — > sg total rate of even just a few 
percent as the significantly harder kaons near Pk = p™ ax that result from the penguin decays 
would dominate over the softer kaons from the ordinary cascade B — > D(+X) — > K(+X'). 

The exclusive B — > Kir decay mode is sensitive to the chromoelectric penguin term. 
Present data and standard model estimates already would appear to exclude enhancement 
of this operator by more than a factor of two or so, although potential uncertainties exist 
concerning the evolution of such "off-shell" operators that may yet modify these estimates. 
Even barring such complications, these conclusions while suggestive nevertheless do not 
preclude a significantly large enhancement in the total, inclusive production of kaons as 
the chromoelectric operator vanishes for on-shell gluons and hence does not contribute to 
the perturbatively calculated, inclusive b —> sg decay rate (obtained using chromomagnetic 
penguin transitions). 

Given these initial results and in view of the fact that enhanced b —>■ sg decays are indi- 
cations of new physics, the simple first step of looking for it via the inclusive kaon spectrum, 
complemented by better data on the exclusive two-body decays, seems an undoubtedly 



1 We are indebted to H. J. Lu for this observation. 
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worthwhile enterprise. 
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FIGURES 

FIG. 1. The expected distribution of kaons (normalized to one) with momenta coming from 
the cascade decay Eq. (||) of the B meson. 

FIG. 2. The expected distribution of kaons if 20 % of the decays arose via gluonic penguins 
and the s — > K fragmentation function was SU(3) symmetric with u — > ir + . In dots, the tail of 
Fig. (1) has also been included for comparison. 

FIG. 3. Same as Fig. (2) except a Heavy Quark fragmentation function has been used for 
s —* K. The various plots, from top to bottom, are for e s = 1,4, 10 respectively, in Eq. (pc|). 

FIG. 4. The diagram with largest imaginary phase contributing to B — > Kir in a perturbative 
QCD analysis. The cross indicates the heavy quark propagator that can go on-shell. The square 
represents gluonic penguin operators from H e ff. 
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Figure 4 



